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It  is  proven  in  this  note,  that  if  weighted  measures  approximating  the 
first  boundary  value  problem  for  the  heat  conduction  and  wave  equations  are 
stable  with  respect  to  the  initial  conditions,  the  dispersion  method  which  is 
used  to  find  the  solution  on  the  upper  layer  is  stable. 

1.  We  will  consider  a  weighted  method,  which  approximates  the  first  boun¬ 
dary  value  problem  for  a  one  dimensional  heat  conduction  equation: 

\h  +  A(ay>+>  +  (1  —  a)y>)  =--  ip,  *  e  (0, 1),  t  >  0, 
jr(0,  o  =  v.(J),  y(i,  t)  V(X,  0)  —  u,(x),  (1) 


and  a  symmetric  method  for  the  wave  equation: 

V7»  + A  +  U  — +  oyi-1)  =  <p,  *e(0, 1),  t>t, 

y(0,0  =  Vi(<).  i/(l ,«)=v2(<),  v(*,0)  =  ut(z),  y(x,t)  ■=  u,(x). 


Here 

Ay  =  — 

XX 


p  . 


v7,- 


y’*1  —  2y>  + 


h  and  t  are  steps  in  space  and  time.  The  stability  conditions  for  methods  (1), 
(2)  have,  because  of  the  initial  conditions,  the  form  (see,  for  example,  [1], 
pp.  567,  582-583) 

i  A* 

<J  is _ - 

2*  4t»  '  (3) 


Here  A  —  l,  2  for  methods  (1),  (2)  respectively.  The  problem  of  finding 

solutions  to  (1),  (2)  on  the  upper  layer  reduces  to  solving  the  first  boundary 
value  problem  for  an  ordinary  three  point  difference  equation 

-y-xx  +  ^y  =  /’  ,6(0,1)l  (4) 

y  (0)  =  vi,  y(l)  =  vj. 
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Efficient  methods  for  multi-dimensional  nonstationary  equations  and  iterative 
methods  used  to  solve  stationary  equations  also  lead  to  a  problem  of  type  (4). 

An  efficient  method  for  solving  problem  (4)  is  the  dispersion  method  {see 
[2],  pp.  281-293),  whose  formulas,  in  the  given  case,  have  the  form 


ft  "*  Oi+iVi+i  + 

fK  —  Vi, 

i  -  o,  1, . 

2  +  hi/axk  —  at  ' 

on  —  0, 

<  —  1.2,. 

..,N-  1 

fit+t  + 

Pi  *■  Vi, 

<-1,2... 

(5) 


The  dispersion  method  is  stable  if  the  condition  1.  is  satisfied. 

However,  a  well-known  sufficient  condition  for  the  stability  of  the  dispersion 
(see  [3],  pp.  386-394)  o>0 

(6) 

narrows  down  the  initial  class  of  methods  for  which  only  condition  (3)  holds. 
Thus  an  important  class  of  methods  with  a  higher  order  of  precision  with 

c  a  -  0.5  —  /  12t  £0r  method  ( 1 ) ,  and  <2t*  (<j.  2*  0.25  -V  /  Or*) 

for  method  (2),  which  are  absolutely  stable  with  respect  to  the  initial  condi¬ 
tions,  may  not  satisfy  condition  (6).  More  general  than  (6),  the  stability 
condition  of  the  dispersion  method  is  given  by 

Theorem  1.  If  methods  (1),  (2)  are  stable  with  respect  to  the  initial 
conditions ,  then  the  dispersion  method  for  finding  the  solution  on  the  upper 
layer  is  stable. 


Before  we  prove  Theorem  1.,  we  will  consider  some  lemmas.  We  will  say 
that  the  dispersion  method  can  be  applied  to  solve  the  three  point  difference 
equation  if  for  the  dispersion  coefficients  the  conditions  |a,l  s£  Ut.  |M  *£  Af,,hold 
where  Mu  Af»  are  some  constants. 


We  will  consider  the  problem 

Ay  =*  —  (ay-),  -f  d(z)y  =  /(x),  *  e=  (0,  1), 

y(0)  >=»  vi,  y(S)«=v,,  o(x)>et>0. 


(7) 


Or,  the  set  S  of  the  grid  functions  v,  which  vanish  at  the  end  points  of  the 
interval  [0,1],  the  tri-diagonal  Jacobi  matrix 


'N  l 


corresponds  to  the  operator  A 


which  holds  for  all  k  =  0,  1,  2, 


Ci  —  At 

0  .  . 

.  0 

0 

0 

—  At  Ct 

—  At .  • 

•  0 

0 

0 

0  —A, 

c,.  . 

.  0 

0 

0 

0  0 

0  .  . 

■  ~  AN-t 

CN-t 

~AN-1 

« 

0  0 

0  .  . 

.  0 

AN- 1 

CH-X 

here  A,  *= 

a  1 1  h\  C, 

“  1  +  di. 

It  is  easily 

related  by  the  recurrence' 

relation 

—  A*+i|A*-i 

. 

(8) 

>  M  1  I 

we  set 

H 

II 

1 

tR-l 

]  =  o. 
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Lemma  1.  The  solution  to  problem 


A»  —  0, 

«’(0)-0,  »(!)- 


*  ^  (0,  l), 

jr 


(9) 


is  positive  if  A  is  a  positive  definite  operator  on  S,  and  has  alternating 
signs  if  A  is  negative  definite. 


Since  A  has  fixed  sign,  the  solution  to  problem  (9)  exists  and  is  unique. 
Using  (8) ,  we  find  it  can  be  written  in  the  form 

« 

D<’m  *  n  I**-*  i- 

«— i 

By  the  well-known  Sylvester  criterion  (see  [4],  pp.  276-279),  the  positive 
definiteness  of  A  implies,  that  all  |rt»|>0,  and  negative  definiteness 
implies  that  the  sequence  has  alternating  signs.  Since  all  A  are 

strictly  positive  the  lemma  is  proved. 


Theorem  2.  The  dispersion  method  is  applicable  to  problem  (7),  if  the 
operator  A  has  fixed  sign  on  the  set  S.  Proof:  The  dispersion  formulas  for 
problem  (7)  have  the  form 


in 
<*(  +  ! 


“  +  Pm, 

■4,+, 

=  C,  —  4, a, ' 


Vi r— vfe  1  —  0,1 . N—i; 

0,  1-1,2,...,*-!; 


P<-M 


4,0,  +  f, 
4,+i 


«i+i. 


P*  —  V|, 


1-  1,2,  ...,*-!. 


(10) 


Therefore,  the  dispersion  method  is  applicable  if  C,  —  4,a,  **  o  for  !<<«;*_! 

We  introduce  the  grid- function  v^  as  follows. 

i>i  —  a<+j(r<4.|,  <  — 0,  1,  ...,*  — 2, 

H 

**  -  *JJ4.-*|**_,|. 

•  Ml 

Using  formulas  (10),  we  see  that  v.  satisfies  problem  (9).  Assume  that 
the  conditions  of  the  theorem  are  satisfied;  then  by  Lemma  1  the  solution  of 
problem  (9)  is  bounded  and  does  not  vanish  at  any  point  of  the  grid.  Therefore 
the  dispersion  coefficient  a.  together  with  the  6.  are  bounded.  It  is 
easily  seen  that  the  sign  or  is  determined  by  the  fcixed  sign  property  of  A. 
This  proves  the  theorem. 

We  shall  now  find  sufficient  conditions  for  the  fixed  sign  property  of  the 
operator  A  on  S. 
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Lemma  2.  If  the  condition  d(*)  >  -$Cl  +  6,,  (11) 

is  satisfied  for  problem  (7)  then  A  5*  4iE,  and  if  condition  12  is  satisfied 

2 

»'(*)<  - pT«(*)+“(*  +  *))-«».  (12) 

then  A  <  -S,E,  where  E  is  the  identity  operator  and  fi».  are  positive 

constants. 

In  fact,  let  y(i)sS.  Using  the  simple  inequality 

and  the  partial  summation  formula,  we  obtain,  when  (11)  is  satisfied 

N- 

(Ay,  yj  =  2  y,Ay,fc  =  (“y ! ,  i  1  +  (4y .  y)  >  ((&>  +  <0  v.  v.)  >  fii  (y,  v). 


i .  e .  A  >  fiifi. 

Further  using  the  Cauchy- Bunyakovski  inequality  and  the  arithmetic- geo¬ 
metric  mean  inequality,  we  obtain 

W_(  N-l 

2  “il/illi-i  |  (at  +  ai+,  )yi*.  (i; 


Suppose  now  that  (1?)  is  satisfied.  Taking  into  consideration  (13),  we  obtain 

.v  -v-i  .. 

(Ay,  y)  =  2  aiyl  ih  +  {dy'  V)=  2  !><“.  +  a«+l)  +  d i  I  y'h ~~ 

i-t  i-1 

-  2  *tviy<-iA  <  2  [-p-  (‘i  +  •«+»)  +  di  J  y'h  <  -  *• (y*  y>* 


i.e.  LM»,y)  < — #»(». »)•  This  proves  the  lemma. 

Mote.  If  We  denote  by  and  the  minimum  and  maximum  eigenvalues 

of  the  operator  A»y then  it  is  clear  that  in  lemma  2  we  can  require 
instead  of  (11)  and  (12)  that  the  inequalities,  <*(*)> —ki  +  61 ,  d(*)  =5  s£  —  XA_i  —  5., 
be  satisfied.  Proceeding  analogously  as  in  (5),  we  can  prove  ih  this  case  that 
v—  ’  Oh ).  where  y  and  y*  are  the  exact  and  approxi  ;ate  solutions 

to  problem  (7)  respectively^ is  the  relative  rounding  error  resulting  from  arith¬ 
metic  operations. 

We  will  now  determine  the  conditions  for  which  the  dispersion  method  will 
be  stable  for  problem  (7). 


-4- 


Lemma  3.  If  the  condition 

<<(*)  2*  0, 

is  satisfied  for  problem  (7)  then  0^a(<i;. 
satisfied 

2 

d(  j)S* - -(“(*)+  flt*  +  *))• 


(14) 

if  condition  15  is 

(15) 


then  — i<o(^0,  and  the  dispersion  method  is  stable. 


In  fact,  condition  (14)  implies  that  Ci  ^  At  +  Al+I,  and  condition 

(15)  implies  that  Ci  <  —  (A,  -f  « -n ) .  From  here,  using  ai  =  0  and  the 

estimates 


*l»  +  l 


C.  -AtOi  1+(A,/A,+1)(1  -a,) 

Ai+i  1 


a,  , i  ^ 


A*U  T  ut)-t"  Ai+t 


1  +  (Ai/Ai+t)  (1  +  a<) 


when  (14)  is  satisfied  and 


when  (15)  is  satisfied. 


we  establish  the  validity  of  the  lemma  by  induction. 


3.  We  will  now  prove  Theorem  3.  For  problem  (4)  the  conditions  of  lemm 
3  relating  to  the  stability  of  the  dispersion  method  take  on  the  form 


1  4 

—  <  —  - — . 
or*  A* 


(14') 

(15') 


Assume  that  conditions  (3)  are  satisfied,  i.e.  the  initial  methods  are  stable 
with  respect  to  the  initial  conditions.  If,  in  addition,  17  >  °.  then  (14') 
is  satisfied  and  the  dispersion  method  is  stable.  If  it  turns  out  that  a <  0, 
then  from  (3)  we  have 


4t*  1  1  4  1  4 

- <  1  B  -  <  - -  - - -  - -  - . 

2  */.»  or*  t*(1/2‘-A74t*)  A*  1  —  4t*/2*A!  A* 


In  this  case  (15')  is  satisfied  and  the  dispersion  method  is  stable.  This 
proves  the  theorem. 

It  is  the  author's  pleasure  to  express  his  gratitude  to  A. A.  Samarskom, 
who  drew  his  attention  to  the  problem  studied  above. 

Submitted  for  publication  1/27/69 
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